Polarized currents in Coulomb blockade and Kondo regimes without magnetic fields 
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We present studies of the Coulomb blockade and Kondo regimes of transport through a quantum 
dot connected to current leads through spin-polanzing quantum point contacts (QPCs). This struc- 
ture, arising from the effect of lateral spin-orbit fields defining the QPCs, results in spin-polarized 
currents even in the absence of external magnetic fields and greatly affects the correlations in the 
dot. Using equation-of- motion and numerical renormalization group calculations we obtain the con- 
ductance and spin polarization for this system under different parameter regimes. We find that the 
system exhibits spin-polarized conductance in both the Coulomb blockade and Kondo regimes, all 
m the absence of applied magnetic fields. We analyze the role that the spin-dependent tunneling 
amplitudes of the QPC play in determining the charge and net magnetic moment in the dot. These 
effects, controllable by lateral gate voltages, may provide an alternative approach for exploring 
Kondo correlations, as well as possible spin devices. 

PACS numbers: 72.15.Qm, 72.25.-b, 72.10.-d, 73.23.Hk 
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I. INTRODUCTION 

Electronic transport in semiconducting nanostructures 
is studied both as it has great potential apphcations 
in spintronicsji and because of its exquisite control of 
parameters, which ahow insightful probes into funda- 
mental physical phenomena. Electrons in such systems 
experience externally controlled confining environments 
that result in strong Coulomb interactions with other 
electrons. As such, quantum dot (QD) structures pro- 
vide well-characterized and defined systems for studying 
quantum many-body physics. QDs may also allow the re- 
ahzation of solid state quantum computation devices as 
well as spintronic semiconductor devices with unprece- 
dented functionalities!^ 

Manipulation of spin-polarized current sources is cru- 
cial in spintronics. This typically requires efficient spin 
injection into conventional semiconductors. The difficul- 
ties with spin injection from ferromagnetic metal leads 
has stimulated extensive efforts to produce spin polar- 
ized currents out of unpolarized sourcesi^ In this con- 
text, the Rashba spin-orbit (RSO) coupling mechanism 
provides a basis for possible device applications)^ This 
coupling, arising from interfacial structure asymmetries, 
depends on the materials used as well as on the confine- 
ment geometry of the structuresi^i^ Most interestingly, 
the Rashba effect allows external tunability, which has 
been studied experimentally in QDs^ and quantum point 
contacts.i^ 

In this paper we study the electronic transport through 
a quantum dot connected to polarizing quantum point 
contacts (QPCs) in both the Coulomb blockade (CB) and 
Kondo regimes. Due to strong spin-orbit interactions;^ 
QPCs can exhibit spin-dependent hybridization of the 
QD states with the leads, without applied magnetic fields, 
opening the possibility for generating spin-polarized 
transport in an all-electrical setup. These effects are 



controllable by lateral gate voltages applied on QPCs, 
resulting in spatially asymmetric structures, as in re- 
cent experiments.— Using the equation-of-motion tech- 
nique and numerical renormalization group (NRG) calcu- 
lations we obtain the electronic Green's function, conduc- 
tance and spin polarization in different system regimes. 
Our results show that both the CB and Kondo regimes 
exhibit non-zero spin-polarized conductance in this sys- 
tem. We analyze how the spin-dependent hybridization 
of the QPC modifies the charge accumulation in the dot, 
as well as the density of states (spectral functions) of the 
system. Interestingly, we find that the polarizing QPCs 
produce spin polarization and split DOS in the Kondo 
regime akin to that reported for current injection from 
ferromagnetic leads)^ although here it occurs for unpo- 
larized reservoirs. Our theoretical studies suggest that 
these effects could be accessible in experiments and re- 
sult in future spintronic devices. 



The paper is organized as follows. First, we present 
the main features of the current polarization in quantum 
point contact systems with lateral spin-orbit, obtained 
by scattering matrix methods in section |lll The next 
section describes the model of a quantum dot connected 
to current leads via polarizing QPCs. In section ITlIAl we 
use the equation-of-motion (EOM) technique to obtain 
the dot Green's function and discuss numerical results 
in both the CB and Kondo regimes. In section IIII Bl 
we revisit the problem using the numerical renormaliza- 
tion group approach, which provides the most reliable 
information in the Kondo regime for realistic system pa- 
rameters. The paper concludes with final remarks and 
summary in section HVl 



2 



II. POLARIZING QPC 



The effective electric field in the z direction that cre- 
ates a two-dimensional electron gas (2DEG) confined to 
the x-y plane results in the well-known Rashba spin orbit 
interaction,i^ 



^so 



{o-xPy - a-yPx) ■ 



(1) 



Here, and (Jy are Pauli matrices, and a is Rashba spin 
orbit coupling constant, which is proportional to the field 
and is therefore material and structure dependent. Elec- 
trons on this 2DEG entering a quantum dot, pass through 
QPCs defined via a confining potential V{x, y) which can 
be thought as made of two parts: U{y) + Vb{x, y), where 
U{y) defines a hard- wall potential of width W (related 
in the experiment to the side-wall etching defining the 
structure) outlining the overall channel structure, while 
Vb{x,y) is the potential generating the QPC barrier, and 
effectively modulated by the side gate potentials of the 
structurci^ We model such barrier bj*ii 



1, 



Vb{x, y) = ^K, [ 1 + cos — ) + ^mw^y^e(y) (2) 
with y = y ~ys, and 
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FIG. 1: (color online) Total and spin-dependent conductances 
for an asymmetric QPC as function of Rashba spin-orbit cou- 
pling a, obtained from a scattering matrix approach.— The 
interplay of vertical and lateral SO effects may result in large 
asymmetry for the up and down spin components for realis- 
tic structure parameters, even in the tunneling regime shown 
here (C <C 1). 
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where Q{x) is the step function, m is the effective mass 
of the electrons, Lx is the unit length of the structure in 
the X direction (along the current direction) and uj is the 
confinement potential frequency. Notice that this poten- 
tial form is asymmetric in the y direction, to rcficct an 
essential ingredient in the experiments: the QPC poten- 
tial must lack y-reflection symmetry in order to generate 
the polarizing effect along the z-directionj^ir— In fact, 
the fields forming Vb generate a spin-orbit coupling given 



^so = -f • (^ X P), 



(4) 



where /? is material-specific. Notice that VVJ, lies in the 
x-y plane, so that the barrier fields induce a lateral spin- 
orbit coupling. The total Hamiltonian of the QPC will 
then be given by. 



H 



p2 



p2 
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2m 



H§o + Vix,y) + V, 



so- 



(5) 



Using a scattering-matrix formalism to study the spin- 
dependent electron transport in this QPC,^^^ we find 
that a net spin-polarized conductance is produced only 
for y-asymmetric potentials.— We can calculate the con- 
ductance of the structure, assuming that the SO coupling 
a and V{x,y) are zero at the source and drain 2DEG 
reservoirs, while both of these terms in the Hamiltonian 



are turned on in the QPC region. Typical structure pa- 
rameters in experiments can be cast in terms of char- 
acteristic length and energy scales, Lq = 32.5 nm, and 
Eo = 3.12 meV, with ao = EqLq = 1.0 x IQ-^^ cVni, 
a typical value of spin-orbit coupling. Using W = ^ 
2Lo and Wi = O.GLq as width/length of the confining po- 
tential, with w = 6 X lO^^s-i, and ^ = 0.97 x IQ-^^m'^, 
gives results as shown in Fig. [TJ This figure shows spin- 
dependent conductances as functions of Rashba coupling 
a, for a potential barrier which is near its conduction 
onset (or "pinch-off', as controlled by the value of Vg); 
arrows indicate the results of spin up and down conduc- 
tances. We see that conductances and Gj, can be very 
different from each other, even in the tunneling regime 
(where each ^ 1) in which the QPC would oper- 
ate to create a quantum dot in the 2DEG. We stress 
that for these realistic values of structure parameters, 
one obtains non-zero spin polarization even when there 
is no external magnetic field and the injection is unpo- 
larized. This interesting result can be understood from 
anticrossing features in the subband energy structure in 
the channel region defining the QPCi^i^ The spin mixings 
and avoided crossings generate spin rotation as electrons 
pass through the narrow constriction of the QPC, and 
can generate large values of the ratio Gf/Gi, even in the 
tunneling regime. Two of these QPCs can then be used 
to define the QD and result in interesting charging and 
conductance regimes, as we will see below. 
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III. QUANTUM DOT WITH POLARIZING 
QPCS 



equations of motion and numerical rcnormalization group 
formalisms. 



In order to address the transport through a quantum 
dot formed with polarizing QPCs, we consider the single 
impurity Anderson model given by the following Hamil- 
tonian: 



Ikcr cr 

+ H t'^i^la^'ika + 4k„Cda), (6) 
Ika- 

where clj^„{cda) is the creation (annihilation) operator of 
an electron of spin a in the dot. The quantities £ik, £d are 
the energies of the electrons in the conduction band 
channels {I = L,R) and the single local energy level in 
the dot, respectively. U is the Coulomb repulsion be- 
tween electrons occupying the QD with Uda — c^^^Cda, 
while represents the lead-QD hybridization occurring 
via tunneling through the QPC, and which is assumed 
to be fc- independent. The density of states for con- 
duction electrons in each lead is taken to be constant, 
p^{e) = pji{e) = p = {l/2D)e{D - \e\), where D is the 
conduction band halfwidtli (hereafter taken as our energy 
unity) . 

The theoretical description of such quantum dot sys- 
tem, especially in the strong correlations regime, has 
been greatly developed over the years.— Techniques of 
note include quantum Monte-Carlo,— equations of mo- 
tion for the Green's functionSfii and the numerical rcnor- 
malization group approachji^ In what follows, we explore 
the role that SO interactions play on the Coulomb block- 
ade and Kondo regimes of transport of the QD, utilizing 

I 



Equation of motion approach and numerical 
results 



To calculate the charge and conductance of the sys- 
tem we calculate Green's functions (GFs), which allow 
us to take into account the correlations induced by the 
Coulomb interaction in the QD. The retarded double- 
time Green's functions are defined as {h = 1)^ 



i{{A-B))^ 



[A(r),i?(0)] + )e(r)e- 



(7) 



where A and B are generic fermionic operators, [A, 
indicates their anticommutator and (• • • ) indicates the 
thermodynamic average for T > 0, or the ground state 
expectation value for T = 0. The GF can be obtained 
using equation of motion (EOM) techniques, so that 

u^{{A-B))^{[A,B] + ) + {{[A,H]-B)), (8) 
where [A, B] represents a commutator. Iteration of this 
formula generates a hierarchy of expressions, starting 
with the local one-particle GF as 



w-Ed- V 2_ ((cd<,;4^)) = l+U{{cd^nda.c\^)), 

\ k '^-^k) 

(9) 

where t„ = v^io- and a = —a. The new (higher or- 
der) GF on the right hand side of Eq. ([9]) can also be 
determined from ([5]), giving 



{Uj-Ed- U){{cdanda;clj) = (uda) + ta {^{{ckarida; c\^)) - {{ckac\gCda; 4^)) + ((4<f CdcrCdai c\^))) ■ (10) 



Coulomb blockade regime 



Although the EOM in Eq. is exact, a solution 

of the impurity GF requires a procedure to truncate 
and/or decouple the higher order terms appearing on 
the right hand side of pU]) . A solution that captures 
the Coulomb blockade physics is given by the Hubbard-I 



approximation 



.20 



{{Ckanda]c\^)) 
{{Ckac\gCda;c\^)) 
{{cla(^daCda]c\^)) 



{nda){{cka\c\^)) 

{c^ksc\g){{cka]C^da)) 

(cLcdff)((cfc^;cL»' 



(11) 



which allows one to write 



Gd^(w) = {{Cda]c}^„))u = 



da I 



1 - 



(12) 



-ea-U 



is the local GF in the 



where G^^(w) 

"atomic" approximation (the exact result for = 0), 
and g{u}) = X]fc(^ ~ ^k)~^ is the non-interacting GF 
of the leads. The DOS of the system (proportional to 
the imaginary part of Gda) contains two Hubbard peaks 
of width proportional to Fg- = nt'^/D, resulting in the 
broadening of the poles of G°^. The spectral weights 
of these peaks are controlled by the dot level occupancy 
with opposite spin, and caused by the Coulomb inter- 
action in the dot. Notice that the SO-induccd polar- 
ization of the QPC results in different peak widths for 
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the different spins. The Hubbard-I approximation pT|) 
is known to be valid for a large U /T ratio, when the 
Hubbard subbands are well separated in energy scale.— 
It is the simplest scheme which describes correlated elec- 
trons, although, since it ignores the Kondo effect, it is a 
reasonable description only at temperatures higher than 
the Kondo scale (T Tr— see next section). 
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The occupancies of spin t and J, are calculated self- 
consistently from the equation 



{nda) 



TT 



(13) 



where /(w) is the Fermi function. It is clear that when 
the QPCs arc not polarizing, — t^, the occupancy 
curves for spin-up and down coincide and a plateau of 
width ^ U appears when the QD level moves below the 
Fermi level [sd + U >^ Ep ^ £d)- This situation changes 
when the QPCs are polarized (see Fig.[2|), as the different 
ta result in F-i- 7^ F^, which in turn produce different {ridi;) 
and (iidi), especially on both sides of the plateau. 

The zero-bias conductance is calculated using a Lan- 
dauer formula generalized for interacting systems^! 



Gn 



de^^lm[Gdaie)] 
oe 



for symmetric coupling of the leads for each spin, 
can also calculate the polarization factor 



(14) 



One 



FIG. 2: (color online) Occupancies for spin-up {rid^}, spin- 
down {udi) and total spin (ud^) + (nai) vs. Ed at zero temper- 
ature. Parameters used are Tf = 0.06, Fj, = 0.03, U = 0.3, 
with D = 1. Notice asymmetry in (jido-) on each side of the 
plateau. 




(15) 



which gives a measure of current polarization in the sys- 
tem. Figure [3] shows the spin-dependent conductances 
and polarization for the system in the Hubbard-I ap- 
proximation. As (n^t) 7^ ('^di) in general (except at the 
particle- hole symmetry point, Ed = — C//2), the conduc- 
tance per spin Go- are also different. Notice that the 
spin-dependent conductance peaks arc very asymmetric 
and non-Lorentzian, due to the peculiar behavior of the 
occupancies and their different up and down-spin cou- 
plings. As we consider here the case Ff > F^, one clearly 
sees that generally G^ > G^ over the entire range of Ed 
values. As a consequence, there is a net up-spin polar- 
ization (~ 60%) and conductance around the resonant 
peaks, with the latter reaching ~ 0.3(e^//i). 



2. Kondo regime 



FIG. 3: (color online) Spin-dependent conductance and po- 
larization as function of Ed at zero temperature. The other 
parameters are as in Fig. [21 



To study the low-temperature behavior of the system 
within the EOM we need to consider higher order GFs 
in Eq. ([T0|) . Using Lacroix's approach,— one can obtain 
relations for the three GFs on the right side of ([TU]) as: 
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k 



(16) 



and 



{Ul - Eda + Sda - eka){{ck9cli^Cda;cl^)) = ( [4^^^^^^ ; 4^] + ) + {{[Cd^CkaCda] H]; clj) 

= (cLcte) +ia{{ndaCda;clj) + Y^hta {(4, ^CkgCda; cl^)) 

k' 

+i'y{{CdgCk9Ck'a;cl„))], 

(W - Eda ~ Sds + Eks - U){{clgCdaCda\ c\^)) = ( ^da Cda ; 4^] + ) + ( ( [4a ; i?] ; 4<t)) 

fe' 

-is{{cl-CdaCk's\c\^))]- 

I 



(17) 



(18) 



Following the decoupling procedure in Ref. [13, each GF resulting in an equation for the dot GF given by 
of the type {{A*BC, D*)) is replaced by 

« {A*B){{C,D*)) - {A*C){{B,D*)), (19) 



Gda(w) = [U{uj)-{nds)-Bs{Lo)-Bs{LOi)] 

X {U{lo)[lO - Sda - + \Bs{i^) + Ba(^l)] Sa(c<j) - As(c<j) + ^^(cJi)}-! , 



(20) 



where I]^(a;) = Efe l^^^l V(a; - Efca), ^ [C/ - a; + e^^ - I]„(cj) + ^^(a;) - (wi)]/t/ and cji = -u; + e^^ + J7. The 

functions Acr{uj) and Bcr{uj) are given by 



■ S^o/^-E^ E:(c.')-Sa(c.) 



(jj — uj' — iS 



uj — uj' + iS 



{i + GdA^')^A^')) 



-(i+GL(^')s:(c.')) 



(21) 



and have to be calculated self-consistently. In the limit 
[/ -> oo, the GF ((20)) acquires a simpler form, 

- e<icr - (1 - -DCT(ti^))Ecr(CL;) - Aa(UJ) 

In the following, we solve numerically for the spec- 
tral function (DOS) in the Kondo regime by the self- 
consistent iteration of Eqs. ([T^ and (|22p . Figured shows 
the DOS vs. u! at low temperature (T = 10~^) for both 
spin orientations (and total), when the QD electron level 
Ed is taken to be at —0.136. The three curves exhibit 
Kondo resonance peaks near the Fermi level (w ~ 0), in 
addition to a much broader peak at w ~ e^jji Since the 
hybridization between leads and QD are spin dependent 



in this case, the DOS clearly splits into two different com- 
ponents for spin up and down. The DOS for spin down 
is shifted upwards in energy with respect to the spin up 
component, resulting in a lower occupancy for the down 
spin. Figure [5] shows indeed the occupancy curves for 
(7icit), (ndi) and total (n) vs. Ed at a given temperature 
(T = 10^^). Qualitatively similar to the results in the 
Coulomb blockade regime, the occupancy curves show 
that {udf) > {ndi) for to > —0.05, while the relation is 
reversed for smaller uj values, reflecting the asymmetry 
introduced by the polarizing QPCs, through r-|- and T^. 
Notice also that the plateaus at 0.5 are not as well defined 
here, due to the enhanced spin and charge fluctuations 
in the Kondo regime. The spin-orbit effect can be under- 
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FIG. 4: (color online) Density of states in QD as a function 
of (xJ in the Kondo regime. Parameters used are T = 10~^, 
Ed = -0.136, Tt = 0.06 and r,^ = 0.03 {D = 1). Sharp 
feature near the Fermi level {u — 0) is the signature of the 
Kondo screening, although asymmetric here for the different 
spin species. 



stood qualitatively as arising from a shift in the dot level 
(as well as depending on the occupancy factors): since 
the effective level position of the electron with spin cr is 
given by Sda ~ Ed + ^^^^(a;), where oc is the 

self-energy, and thus naturally causes the spin-dependent 
occupancy seen in the figure. 
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FIG. 5: (color online) Occupation in the QD as function of 
Ed. Parameters used as in Fig. |?1 

Figure IH^a) shows the spin-dependent conductance 
curves vs. e^. Several features are noteworthy. As 
Ed changes, the spin-dependent conductances exhibit the 



FIG. 6: (color online) (a) Conductance, and (b) polarization 
and difference between spin up and down conductance, as 
functions of ed- Parameters as in Fig. [H 



anticipated peaks at low temperature, with spin up con- 
ductance dominating (naturally, as > F^). Figure 
ini^b) shows the difference between spin up and down con- 
ductances as well as the spin polarization vs. Ed- In this 
regime, the difference between the conductance for the 
two spin orientations reaches ~ 0.4(e^//i). Correspond- 
ingly, the net spin polarization reaches rj ~ 60%. 

Let us now analyze in more detail the effect of temper- 
ature on the conductance of the system and especially 
its drop for Sd ^ 0. The conductance curves in Fig. 
inja) are highly asymmetric about the Fermi energy and 
vanish rapidly away from it. This vanishing for very neg- 
ative values of Ed is due to the the Kondo temperature, 
Tk, becoming smaller than the temperature of the sys- 
tem. Figure [7] presents the total zero bias conductance 
as function of for several values of temperature. As T 
is lowered, the total conductance increases for a given Sd, 
and the width of the conductance peak increases in Sd- 
This explicitly reflects the existence of the Kondo reso- 
nant peak in the spectral function, and how the system 
will reach the unitary limit of conductance for T = for 
Ed well below the Fermi level. 

Our discussions above arc applied to the case of infinite 
U, where the EOM method gives qualitatively accurate 
results in the Kondo regime. Extension of this approach 
for finite U is known to be problematic, including its 
failure to exhibit a Kondo resonance at the particle-hole 
symmetry point (e^ = —{7/2).— To carry out our study 
in the finite U case, we use instead the essentially exact 
numerical rcnormalization group approach, as we discuss 
in the following section. 
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FIG. 7: (color online) Total conductance as function of for 
several temperatures T . Notice that Tk is strongly suppressed 
for more negative Ed values, which lowers the conductance at 
a given T . Parameters as in Fig. [l] 



B. NRG results for finite U 




FIG. 8: (color online) Density of states in the QD as function 
of energy obtained from NRG calculations. Polarizing QPCs 
generate effective splitting of the DOS for different spins. Sys- 
tem parameters used are U — 0.5, — 0.06, = 0.03, 
L» = 1, and T = 0. 



For the finitc-C/ case we study the spin polarized con- 
ductance using the standard numerical renormalization 
group approach i^^i^^ In this case, unlike the previous 
infinite-?/ case, the processes involving double occupied 
states are naturally present in the dynamics of the sys- 
tem, and allovif for a reliable description of the low-energy 
behavior. We set U = 0.5 and T = 0. Figure [8] depicts 
the local density of states calculated with NRG for the 
same system parameters as before, F^ = 0.06, Y ^ = 0.03, 
and with e^; = —0.2. This value of corresponds to a sit- 
uation where the system is away from the particle-hole 
symmetric point (= — [7/2), more suitable to compare 
to the previous infinite- ?7 calculation (where the system 
never reaches the symmetry point). Notice that there 
is a strong spin asymmetry in the DOS; on the nega- 
tive side of the w-axis, the DOS for spin up (solid curve) 
presents a peak near —0.2, corresponding to the energy 
of the local bare orbital e^, slightly shifted by the real 
part of the proper self-energy. On the positive side of 
the w-axis, the peak near td~^'^ would result in the suc- 
ceeding CB peak, which appears only for the spin-down 
DOS (dashed curve). In this regime as well, the fact that 
F.f > Tj, favors the spin-up occupancy to the detriment of 
the spin-down occupancy. One also notices the important 
peaks close to the Fermi level, signature of the Kondo ef- 
fect, which are slightly split away and suppressed by a 
seemingly effective magnetic field induced by the spin- 
asymmetric coupling to the leads. This phenomenon is 
akin to the suppression discussed by Martinek et al.^ in 
the context of a QD coupled to ferromagnetic leads. No- 
tice, however, that no external magnetization is present 
in our system and that the polarization is only arising 
from the QPCs and due to the lateral SO interaction. 



Analogously to the infinite-t/ case, the spin-asymmetry 
discussed above induces spin polarized transport in the 
system. In Fig. [3] we show the conductance as a func- 
tion of Ed for the same parameters as Fig. [H] Notice that 
away from the p-h symmetric point [sd = —0.25) the con- 
ductance for spin up is much larger than for spin down, 
resulting in a sizeable polarization (jj ~ 70%), as shown 
by the (green) dotted curve. At the p-h symmetric point 
the conductance for both spins reaches the unitary limit 
and 77 — > 0; this is consistent with the restoration of the 
Kondo state of the system at the p-h point for ferromag- 
netic leads 1^ 



IV. SUMMARY 

In summary, we have investigated the spin-dependent 
transport properties of quantum dot structures with po- 
larizing quantum point contacts. We have shown that as 
QPCs can generate finite spin-polarized currents, due to 
the combination of lateral and perpendicular spin-orbit 
interactions, they also induce current polarization in 
quantum dots made with these QPCs. Using equation-of- 
motion techniques and numerical renormalization group 
calculations, we obtained the electronic Green's function, 
conductance and spin polarization in different parame- 
ter regimes. Our results demonstrate that both in the 
Coulomb blockade and Kondo regimes, the quantum dot 
exhibits non-zero spin-polarized conductance, even when 
the injection is unpolarized and there are no applied mag- 
netic fields. The spin-dependent coupling is shown to give 
rise to nontrivial effects in the density of states of the sin- 
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glc QD, resulting in strong modification of the charge dis- 
tribution in the system. Most importantly, these effects 
are controllable by lateral gate voltages applied to the 
QPCs, and together with the ability to create quantum 
dots, they provide a new approach for exploring spin- 
tronic devices, spin polarized sources and spin filters. 
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